Segregation in noniteracting binary mixture. 



O 

o 

(N 

< 



o 

e 

I 



O 
o 



> 

ON 

cn 

p 

od 

o 
o 



X 



Filip Krzyzewski and Magdalena A. Zaluska-Kotur 
Institute of Physics, Polish Academy of Sciences, Al. Lotnikow 32/46, 02-668 Warsaw, Poland 

Process of stripe formation is analyzed numerically in a binary mixture. The system consists of 
particles of two sizes, without any direct mutual interactions. Overlapping of large particles, sur- 
rounded by a dense system of smaller particles induces indirect entropy driven interactions between 
large particles. Under an influence of an external driving force the system orders and stripes are 
formed. Mean width of stripes grows logarithmically with time, in contrast to a typical power law 
temporal increase observed for driven interacting lattice gas systems. We describe the mechanism 
responsible for this behavior and attribute the logarithmic growth to a random walk of large particles 
in a random potential due to the small ones. 

PACS numbers: 68.35.Fx, 61.20.Ja, 64.60Qb, 5.50.+q 



I. INTRODUCTION 

Binary mixtures subject to an external driving force 
are often found to segregate and form stripes of width 
increasing with time [l|, @j S 0] ■ Such phenomenon has 
important practical applications in chemical or pharma- 
ceutical technology. Much theoretical and experimental 
effort has been devoted to understandingthe main mech- 
anisms responsible for this process [1, [1,13| ■ The simplest 
and best known model of stripe formation under external 
driving force is a simple lattice gas with attractive near- 
est neighbor interparticle interactions. It has been known 
for a long time that such system orders under an influ- 
ence of the external driving force 8, 9, 10]. Recently, this 
model has been studied intensively [ill, [ij, El [IE Ell 
in a context of latest experimental results. We describe 
below the process of stripe formation in a lattice model 
with two types of particles which do not interact via di- 
rect forces. The model is defined in such a way that larger 
particles, which occupy five lattice sites each, can overlap. 
Hence they block smaller number of sites, when they are 
close together, than when they are separated. Presence 
of smaller particles induces effective interactions between 
large particles [l^ and the system orders at high enough 
densities. In the ordering process, due to the presence 
of smaller particles in the spaces between stripes, large 
particles realize random walk in a random potential. A 
jump can occur only when there is enough free space in 
a chosen direction; the latter being a random event. As 
a result kinetics of stripe growth in our model is different 
than typicalpower law temporal growth known in driven 
lattice gas ^ Es'] and a mean stripe width in the binary 
mixture studied here increases logarithmically with time. 



II. STRIPE FORMATION IN THE BINARY 
SYSTEM 



Stripe formation can be easily observed in a simple 
lattice gas model with nearest neighbor attractive inter- 
actions. It has been shown [1, [^, [l3| by means of Monte 
Carlo simulations that such system orders successively 
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FIG. 1: Jump rules for large (at left) and for small (at right) 
particles. Shaded area sites can be shared by overlapping 
large particles. The external bias h affects the jump rates of 
small particles only. 



in stripes under an influence of external driving force. 
When the evolution starts from some random configu- 
ration, stripes arc formed: initially thin they gradually 
become thicker. This process, described and analyzed in 
detail in Refs EH, [U, M: [jJ, [l5|, M occurs in two stages 
- stripe formation and then stripe growth. The mean 
width of stripes grows typically as a power of a simu- 
lation time (where time is measured in the number 
of Monte Carlo steps). Typically two different powers 
X = 1/3 and x = 1/4 are observed. 

Below we present segregation process in a binary mix- 
ture of noninteracting particles. The system consists of 
large objects which can overlap when they occupy neigh- 
boring sites and a large number of smaller objects. Large 
particles occupy area of five sites. The core site in the 
middle is always occupied by a center of one large parti- 
cle only, while the other four sites can be shared by other 
large particles. Small particles occupy one site and do 
not share it with any other particles. When large parti- 
cles overlap they leave more free sites for small particles 
so the number of their possible configurations increases, 
effectively creating forces between large particles. The 
system orders forming stripes, when driven by either pe- 
riodic in time or static external force. The large particles 
move randomly with the same jump probability in each 
direction. At first, we choose randomly one of four di- 
rections. Then occupation of sites in this direction is 
checked. Each large particle moves only if other parti- 
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FIG. 2: Successive stages of separation of large (dark) and 
small (light) particles 
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FIG. 3: Number of free sites in the system as a function of 
time for one sample. System size is (25 x 250) sites populated 
by 500 large and 3972 small particles. The initial free sites 
fraction is 0.01. 



cles (large or small) do not block it. Small particles are 
affected by external, biasing field characterized by a con- 
stant b. Jump rates of small particles in each direction 
depend on b. When the particle moves in the direction 
of the bias the jump rate is p = 1. When it moves in the 
direction perpendicular to it then p ~ b^^. For a jump 
opposite to the field we have p = b^^. In most of our 
calculations we have used 6 = 5. Fig 1. illustrates the 
jump rules for large and small particles. 

There is no direct interaction between particles, but 
they can partially or fully block sites, which then can- 
not be occupied by other particles. Large particles can 
overlap; so they block fewer sites, when they are closer 
together. When the large particles overlap, then the freed 
sites become immediately available to small particles. 
Once they are filled the large ones cannot move away 
from each other. Therefore, a configuration of closely 
packed large particles is more probable than other con- 
figurations. It is shown in Refs IitI, [isl how to map such a 
mixture of large and small particles onto an Ising model 
with nearest-neighbor interactions. Strength of this ef- 
fective interaction depends on the density of the smaller 
particles. The interaction is stronger when the density 
is higher and, eventually, above certain critical density a 
phase transition occurs. This entropic interaction along 
with the biasing field causes a formation of stripes paral- 
lel to the field. The system structure at different evolu- 
tion times is shown in Fig. 2. Time for all presented data 
is measured in a number of Monte Carlo (MC) steps. 

Increasing fraction of free space is a signature of an 
ordering process. Ratio m of the number of free sites to 
their total number as a function of time t is plotted in 
Fig. 3. At t > 10"* the curve shows a step like structure. 
Each step lasts a relatively long period of time (note log- 



arithmic time scale in Fig. 3) during which the number of 
stripes does not change. The steps can be observed when 
there are only few stripes present in the system. A transi- 
tion from one step to the other means that the number of 
stripes is reduced by one. Average number of stripes de- 
creases with time t and their average width l{t) increases 
with time, as seen in Fig. 4. We note that, as in Ref. 
|12| there are two stages of stripe width growth. The first 
stage starts when clusters of large particles coarsen and 
ends at the time when stripes start to become as long as 
the system width and a multistripe order becomes clearly 
visible. The second stage begins with well formed stripes 
and continues through the process of joining them to- 
gether. These two stages can be seen in Fig. 3. The first 
stage of stripe formation corresponds to a line of lower 
slope and the second stage begins when the line bends up 
and becomes steeper. During the first stage the system 
looks like in the first left hand side panel in Fig 2, and 
the second stage is represented in the next four panels. 
Number of unoccupied sites increases with time. 

We carried out simulations for various lattice sizes with 
different width-to-height ratios. Results of all simula- 
tions were averaged over 100 realizations. Large particles 
were distributed randomly at the beginning of every run. 
They could overlap. Then, given large fraction of remain- 
ing free sites was occupied randomly by small particles. 
Then, a biasing field was turned on and the simulation 
started. We have compared results obtained for closed 
systems with fixed particle number with those for open 
systems with fixed chemical potential controlling the par- 
ticle density. Formation of structures has been observed 
at a static bias field and at a field with periodically chang- 
ing orientation. The highest formation rate was noted for 
a constant field, and it decreased with an increasing fre- 
quency of the field variation. The results do not change 
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FIG. 4: Average width of stripes d — kI as a, function of time, 
scaled by parameter k. See text for explanations. 



qualitatively until a frequency of around 1/(5MC steps) 
is reached, above which stripes stop to form at all. Thus, 
most of the presented examples here are calculated for a 
constant in time driving field. 

To find an average width of the stripes, a correlation 
function fc{r) for a given configuration of particles was 
calculated. It is defined as 
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where |r| is the distance between i-th and j-th site along 
the direction perpendicular to the external field; n^, nj = 
0, 1 are occupations of i-th and j-th site, = 1 when site 
i is occupied by one or more large particles and rii — 
when the site is empty. N is the number of large particles. 
Sums are over all sites whose coordinates perpendicular 
to stripes differ by \r\. Average width of the stripes is 
such \r\ for which of the correlation function has the first 
minimum. 

Our simulations show that mean width of stripes grows 
as logarithm of time / = log(t/i||), where i|| is the sys- 
tem width, parallel to stripes — 25 in Fig. 2). The 
character of growth does not depend on the system size. 
This is illustrated in the Fig. 4, where seven data sets, 
are plotted in two groups for the binary mixture and top 
plot is for a one component, interacting driven lattice gas. 
The main panel shows the mean width of stripes I as a 



function of log(t). Three lower data sets, plotted by -|- 
and X , represent results for systems of constant number 
of smaller particles, and four data sets in the middle, plot- 
ted by ■, o, • and A represent results for open systems, 
with fixed chemical potential. For all these plots stripe 
width I was rescaled by a parameter k chosen in such a 
way that each data set lies on one line log[i/(a-L||)] with 
a = 10 for closed and a = 1.4 for open systems. Scal- 
ing parameters for closed systems are k = 0.4, 0.5 and 
0.7 for the ratio r = L\\/L±_ equal, respectively to 2,4 
and 10. System sizes are (50 x 100), (25 x 100), and 
(25 X 100) and number of large particles N = 600, 300 
and 600. For open systems n = 0.833, 1, 1.66 and 1.25, 
respectively. Size of each system is (25 x 100), N — 600 
and 80%, 85%, 90% and 95% sites are occupied, respec- 
tively. The upper set printed by triangles V represents a 
stripe growth for one component interacting driven sys- 
tem at a temperature Q.STc°° , where Tc°° = 3.18J/fci3 
(where J is interaction strength) and jump probabilities: 
p = 1 in the direction of field, p — Q'va. the opposite direc- 
tion, and p = exp{—cJ/T) in the direction perpendicular 
to the field, where c is a number of nearest neighbors. 

We can now compare results for a one-component in- 
teracting driven system with those for a non-interacting 
binary mixture. Stripe formation is observed in both 
cases. It can be seen, however, that time dependence of 
these processes has a different character. In contrast to 
typically noted power law in time stripe width growth 
for one component interacting driven systems, in binary 
mixtures with entropic interactions we typically observe 
slower logarithmic time dependence. Ordering in the bi- 
nary system happens due to the indirect, effective at- 
tractive interaction between particles. Strength of this 
entropic interaction decreases with increasing number of 
free sites [l3]- In the closed system with a fixed parti- 
cle number the available free space expands with time 
(Fig. 3) and the effective interactions weaken. In order 
to check if this weakening has influence onto the charac- 
ter of stripe growth, we studied also open systems with 
varying small particle number and where mean number 
of free sites was controlled by an external potential. As 
illustrated in Fig. 4, the time dependence in open sys- 
tem has the same logarithmic character observed for the 
closed system. We see that the decreasing with time in- 
teraction strength is not the main reason for the type of 
observed temporal width growth. 

Specific aspect of the system studied here is the ex- 
istence of particles of two types. In order to cross an 
interstripe distance, large particle have to find their way 
between densely packed set of small ones. To execute a 
jump, the large particle has to wait until a passing stream 
of small particles creates a hole, large enough to fit in. 
As a waiting time for a jump in such case varies from 
one event to the next, we can treat such process as a 
random walk in a random potential. In the next section 
we show that the logarithmic character of the temporal 
stripe width growth can be explained by such a descrip- 
tion of large particles kinetics. 
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III. MECHANISM OF STRIPE GROWTH 

Stripe growth is an anisotropic process that takes place 
in the driven systems. The main course of growth hap- 
pens along the direction perpendicular to stripes. Ex- 
istence of the second dimension controls relative proba- 
bilities of several mechanisms that compete in the stripe 
growth process. This process in one component system 
has been analyzed and explained in details in Refs [T^ 
and 13. Ref. [ij describes two different competing mech- 
anisms: evaporation/condensation of particles from the 
surface of the stripe and diffusion of particles/holes be- 
tween interfaces. The former one leads to I ~ t-^^^ and 
can be observed at earlier times or for shorter systems, 
whereas the latter leads to the I ^ t^^^ growth and is 
activated at later stages of stripe formation or in longer 
systems. 

Let us consider a one component system with particles 
attracting each other. System orders under influence of 
a static bias field, initially forming many thin stripes. 
When the process continues some of stripes disintegrate 
while the remaining ones become thicker. Stripe extinc- 
tion is a process consisting random actions of a single 
particle: the particle evaporates first from the stripe wall, 
then it walks randomly in an empty space until it read- 
sorbs at the same or the other wall. The process con- 
tinues until one of stripes disintegrates. Decay of one of 
two neighboring stripes is a problem similar to that of 
the gambler ruin. We are not asking, however, a stan- 
dard question about the probability of a ruin. Instead, 
we are interested in the mean time of ultimate decay of 
the first or the one of two neighboring stripes. This time 
is proportional to the mean time of evaporation of one 
particle row across one stripe. Number of particles in 
such row is equal to the width of the stripe and fluc- 
tuates as particles escape from and stick to the domain 
walls. Emergence of a fluctuation of size I means that row 
of such length disappears. Mean time of such an event 
scales as P. Fluctuations occur independently in each 
row, so the time in which the entire stripe disappears is 
proportional to the number of rows in one stripe, Ly and 
to the time tq needed for a single particle to pass the 
distance from one stripe to another. Thus we have 



(2) 



where is time scale parameter. 

Time tq is mean first-passage time of a distance be- 
tween stripes by a particle. In a general case of an in- 
homogeneous potential the first-passage time is given by 
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where s = l{l — p) / p is the distance between stripes, p a 
density of large particles in the system, pj is jump rate 
from site j in the direction pointing from the initial site 



to the final site s, and qj is jump rate in the opposite 
direction. For the one component system we can assume 
that pj = Qj in the Eq. ([3]) and that pj are the same for all 
j — 1, ...s except when j = for a jump originating at a 
site neighboring to the stripe. The rate pq is a probability 
rate for a particle jumping out of other particles. The 
particles attract each other, so this rate is smaller than 
all others: po < pi- The interaction and so the value of 
Po depend on the number of neighboring particles. We 
can write 
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and treat po as an effective rate averaged over many 
jumps. The above formula is correct when the density 
of particles between stripes is low, not higher than one 
free particle per row. If s > pi/po then the second term of 
(HI dominates. In this case, however, the particle density 
is higher, than one particle per row so proper expression 
for the time tq is obtained by dividing by a number of 
particles that reach the wall per time unit. This number 
is proportional to the distance s, so for s > pi/po we get 
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Equation ([5]) is valid if the density of particles between 
stripes is higher than one particle per row but is still quite 
low. For higher densities, however, pair interactions in 
the empty space start to play a role, causing the entire 
process to slow down. Using Eq. ^ we obtain 
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Solution of this equation for To given in Q is 
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(7) 



the exponent of the power law growth of I changes be- 
tween 1/3 and 1/4. When (O is used, we get 



(8) 



i.e. the power law time dependence with a single expo- 
nent a; = 1/3 (l3| . For higher temperatures when particle 
density between stripes becomes higher, various values of 
exponent x, usually smaller than 1/3, are observed. Still, 
I ~ t^/^ is a dominant behavior for wide range of tem- 
peratures and system geometry parameters. 

When the system consists of two different types of par- 
ticles, random walk from one stripe to another is not free. 
Each particle has to wait until there is enough space for 
it jump. We can treat the process of particle motion in a 
dense medium as a random walk in a random potential 
landscape. A jump to the left with rate pi and jump to 
the right with rate Pr are in this approach independent 
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events, occurring according to the same probability dis- 
tribution. Such a model leads to the following expression 
for the mean first-passage time jT§| 
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with being an average over random variable realiza- 
tions. Thus all linear in s terms in expression ^ are for 
large s irrelevant and we get A = log(7)(l — p)/p. Using 
now ^ and ^ we get the following equation 
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Its solution for large I and t can be written as 
l^log{t/Lu) 



(11) 



(12) 



and, indeed such character of the time dependence is ob- 
served in Fig 4 for binary systems. It can be seen in the 
inset of Fig 4. , that power law cannot be fitted to the data 
sets for binary mixtures. The character of stripe growth 
is the same for closed system, where number of free sites 
increases as it is for an open system with constant density 
of small particles controlled by external potential. 



IV. SUMMARY 



We have investigated binary mixture system driven by 
an external force. Particles in this binary system do not 
interact with each other directly but they effectively do 
so via indirect entropy interaction. The system orders 
forming stripes, similarly like in driven single component 
system with attractive forces. 

The existence of two different particle types leads 
to the logarithmic temporal growth of the mean stripe 
width. Such time dependence is slower than the power 
law temporal growth in an interacting one component 
system. In binary systems large particles travel among 
densely packed small particles, which effectively slow 
down their wandering. We have attributed the logarith- 
mic growth process to a random walk of large particles 
in an effectively random potential. 
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